Optical realization of relativistic non-Hermitian quantum mechanics 
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Light propagation in distributed feedback optical structures with gain/loss regions is shown to pro- 
vide an accessible laboratory tool to visualize in optics the spectral properties of the one-dimensional 
Dirac equation with non-Hermitian interactions. Spectral singularities and VT symmetry breaking 
of the Dirac Hamiltonian are shown to correspond to simple observable physical quantities and 
related to well-known physical phenomena like resonance narrowing and laser oscillation. 
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Introduction. Quantum mechanics prescribes that the 
Hamiltonian H of a physical system must be self-adjoint. 
Since the the seminal paper by Bender and Boettcher 
[l|, it was realized that the Hermiticity of H can be 
relaxed, and that a a consistent quantum theory can 
be constructed for a broader class of Hamiltonians 0- 
5]. Among these are parity-time (VT) Hamiltonians, 
which possess a real- valued spectrum below a symmetry- 
breaking point. Non-Hermitian Hamiltonians are also en- 
countered in reduced descriptions of open quantum sys- 
tems, with important applications to atomic, molecular 
and condensed-matter physics @. In such systems, the 
lack of Hermiticity can lead to the appearance of excep- 
tional points and spectral singularities, whose physical 
relevance has been discussed by several authors (see, e.g., 
041). Recently, non-Hermitian extensions of relativis- 
tic wave equations [Iol - [l2| and non-Hermitian quantum 
field theories [HI have attracted an increasing interest as 
well. As some issues in this field are still debated (see, 
e.g., [Hj]), physical realizations of non-Hermitian rela- 
tivistic models remain mostly unexplored. Recently, op- 
tical structures in media with a complex refractive index 
have been proposed to test and visualize non-Hermitian 
features rooted in the non-relativistic Schrodinger equa- 
tion with a complex potential d, H, [3 EH ■ The main 
motivation in the study of such quantum-optical analogs 
is that concepts like exceptional points, spectral singu- 
larities and VT symmetry breaking become measurable 
quantities in an optical experiment p| [9[. This has 
lead to the first experimental visualization of exceptional 
points and VT symmetry breaking in an optical struc- 
ture , E3] • Such results motivate the search for optical 
simulators of non-Hermitian relativistic wave equations, 
which is the aim of this Letter. Here it is shown that light 
propagation in distributed-fcedback (DFB) optical struc- 
tures with gain and/or loss regions, which is at the heart 
of such important devices as DFB semiconductor lasers 
[H, Qj| , can provide a fertile ground to test the spectral 
properties of non-Hermitian Dirac Hamiltonians. Simi- 
larities between light propagation in DFB structures and 
relativistic wave equations were noticed in early studies 
on gap solitons in connection with the massive Thirring 
model of field theory [2(| , however these previous studies 
did not consider non-Hermitian interactions. 
Non-Hermitian Dirac equation and its optical realiza- 



tion. Let us consider the Dirac equation in one spa- 
tial dimension for a two-component spinor wave func- 
tion ip(x,t) — (ipi, '02) T with time-independent vector 
(V) and scalar (S) couplings, which in natural units 
(ft = c=l) reads pj], IS 

idfip = —iad x 4> + j3m(x)ip + V(x)ip = Hip (1) 

where a and /3 are 2x2 Hermitian square matrices sat- 
isfying the relations a 2 = /3 2 = 1 and a/3 + (3ct = 0, 
V(x) is the time-component of a Lorentzian 2- vector po- 
tential, m(x) is the space-dependent effective mass de- 
fined by m(x) = m$ + S(x), and mo is the rest mass 
of the Dirac particle. Among the various representa- 
tions of the Dirac equation, the optical realization of 
Eq.(l) discussed below is at best highlighted by tak- 
ing a = o~ z and j3 = o~ x , where a x = f J J ] and 



are the Pauli matrices. If either V or 
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m assume complex values, the Dirac Hamiltonian H is 
non-Hermitian and its spectrum is, in general, complex- 
valued. For the Dirac equation, parity and time-reversal 
operators can be defined as 12] V = Vq/3 and T = K./3, 
where Vq changes x with —x and /C performs complex- 
conjugation. Hence VTip(x) — ip*(—x). VT invariancc 
of H, i.e. [VT,H] = 0, requires V(-x) = V*(x) and 
m(— x) = m *{x) |12|. If every eigenfunction of a VT- 
invariant Hamiltonian is also an eigenfunction of the VT 
operator, the VT symmetry of H is said to be unbroken 
and its spectrum is real- valued In the following, we 
will consider a real- valued effective mass m(x), whereas 
non-Hcrmiticity is introduced by allowing the vector po- 
tential V(x) to be complex- valued. An optical realization 
of the Dirac equation (1) is provided by propagation of 
light waves in an effective one-dimensional DFB structure 
[ia |. Let n(z) = n - Anh(z) cos(27rz/A + 20 (z)) be the 
effective index grating of the dielectric structure, where 
no is the modal refractive index in absence of the grating, 
An <C no and A are the peak index change and the nom- 
inal period of the grating, respectively, and h(z), 28(z) 
are the normalized amplitude and phase profiles, respec- 
tively, of the index grating. The periodic modulation of 
the refractive index leads to Bragg scattering between 
two counterpropagating waves at frequencies close to the 
Bragg frequency u>b = 7rc/(Ano), where c is the speed of 
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light in vacuum. The linear space-dependent absorption 
coefficient of counterpropagating waves in the structure 
is indicated by ao(x) («o > in lossy regions, «o < in 
gain regions). In a semiconductor DFB structure, gain 
and loss regions could be tailored by a judicious control 
of current injection across the active layer f]~oj ] . Indicat- 
ing by E(z,t) — ip 1 (z,T)exp[-iuj B T + ik B z + iO(z)] + 
tp2(z, t) exp[— iuJBT — iksz — iO(z)] + c.c. the electric field 
propagating in the DFB structure, where ks — 7r/A, the 
envelopes ipi and ^2 of counterpropagating waves satisfy 
coupled-mode equations [19j . After introduction of the 
scaled space and time variables x = z/Z and t = t/T, 
with Z = 2uq A/ (n An) and T — Z/v g , where v g ~ c/uq 
is the group velocity of light at frequency wg, the en- 
velopes ipi and "02 satisfy Eq.(l) with a real mass m and 
a complex-valued vector potential V defined by 

dO 

m(x) = h(x), V{x) = — -i 1 {x), (2) 

where 7(2;) — Zolq(x) is the dimensionless absorption 
coefficient. In the following, we will assume that both 
m{x) and V(x) have a limited support, i.e. m = V = 
for \x\ > L/2. This case typically applies to optical DFB 
structures, in which the grating region is spatially con- 
fined to a finite region of length L. 

Spectral singularities, bound states and resonances of the 
Dirac Hamiltonian. The spectral properties of the non- 
Hermitian Dirac Hamiltonian H, defined by Eqs.(l) and 
(2), can be investigated by standard methods of scatter- 
ing theory. To this aim, let us introduce the functions 
(/>e\x), (j)g\x), tpg (x) and ip^\x), which satisfy the 
equation Htp — Etjj (E is a complex-valued parameter) 

with the asymptotic behavior tf>^ = (1, 0) T exp(iEx), 
= (0, l) T exp(-^x) for x < -L/2, and y§ = 
{1,0) T exp(iEx), ip§ = (0, 1) T exp{-iEx) for x > L/2. 
As the Wronskians W{<j$ = W {ip^ , ip^ } = 

1 do not vanish, {4>e anc ^ W^e^^e} arc two 

sets of linearly- independent solutions to the equation 
(E — H)ip — 0, and therefore there exists a 2 x 2 ma- 
trix M(E), with detM(E) = 1, such that 

4> ( e\x) = M 11 (E)^\x)+M 2 i(E)^\x) 
4 2) (x) = M 12 (E)^\x)+M 22 (E)^ ) (x). (3) 

Physically, M.(E) is the transfer matrix that connects the 
amplitudes of forward- and backward-propagating waves 
from x = —L/2 to x = L/2, which is commonly adopted 
in the optical context (see, for instance, The spec- 

tral transmission and reflection coefficients, for left (I) 
and right (r) incidence, can be expressed in terms of the 
transfer matrix elements in the usual way [s| 
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For a 'PT-invariant Hamiltonian, one has (p^l (x) 
(t> [ E ] *{-x) and (P ( eI(x) = (t> [ E ] *{-x), and thus M^iE) 



(4) 



M*(E*), which implies Mn(E) = M* 22 (E*). 
The spectrum of H , as well as the existence of spec- 
tral singularities arising from the non-Hermiticity of H 
El [HI i can be determined by an inspection of the singu- 
larities and branch cuts of the resolvent G(E) — (E — 
H)~ x , which takes the integral form [2l| G(E)ij){x) = 
J dyQ(x,y; E)tp(y), where Q{x,y;E) is the Green func- 
tion. Its explicit form reads Q(x,y;E) = Q + (x,y; E) for 
lm{E) > and Q(x,y;E) = g^(x,y;E) for Im^) < 0, 
where 



G+(x,y;E) 



M 2 2{E) 



<S>{y-x)^\x)^ )T {y) + 
^(x~y)^\x)^ )T (y)]a x ,(5) 



G-(x,y;E) = 



Mu(E) 



<f(y-x)^\x)^ T (y) + 



(6) 



and $(x) is the Hcavisidc function = for x < 0, 

= 1 for x > 0]. On the basis of Eqs.(5) and (6), the 
following results hold for the spectral properties of H. 

(i) Continuous spectrum. The continuous spectrum of H 
is the entire real energy axis (— 00 < E < 00), where 
G(E) has a branch cut. 

(ii) Point spectrum. The zeros of A4 22 (E) on the 
Im(.E) > half plane, together with the zeros of Adn(E) 

on the Im(E) < half plane, define the point spectrum of 

(2) 

H; at such energies, the function ^ (x) [for Im(E) > 0] 

and <j} ( E } (x) [for lm(E) < 0] are bound states of H . 

(iii) Resonances. The zeros of M 22 (E) [M\x{E)\ on the 
Im(.E) < [Im(-B) > 0] half plane, i.e. the poles of 
the analytic continuation of Q+ [Q-] on the Im(_E) < 
[Im(.E) > 0] half plane, correspond to the resonances [an- 
tiresonances] of the scattering problem. 

(iv) Spectral singularities. A spectral singularity 21] is 
any zero E = Eq on the real axis of either M.n(E) or 
M- 22 {E), around which G{E) is unbounded (for either 
E = Eq + i0 + or E = E a — i0 + ) in spite of the fact that 
E = Eq does not belong to the point spectrum of H . 
Note that, for a 'PT-invariant Hamiltonian, M.n(E) = 
M2 2 {E*) and hence at a spectral singularity both At 11 
and A-^22 vanish, the resolvent is unbounded for E = 
Eq±M + , and the transmission and reflection spectral co- 
efficients diverge according to Eq.(4) |22| . For such a rea- 
son, in Ref. 0] spectral singularities were identified with 
zero-width resonances. However, for a PT-non-invariant 
Hamiltonian, a spectral singularity could arise from the 
vanishing of Aiu but not of A4 2 2, i.e. it might not cor- 
respond to a divergence of the spectral transmission or 
reflection coefficients. Such a case was not noticed in 
Ref.Q and has a different physical meaning: at such a 
spectral singularity the potential becomes reflectionless 
for simultaneous excitation with two waves of same en- 
ergy Eq and amplitudes 1 (for the wave incident from the 
left side) and M.21 (for the wave incident from the right 
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FIG. 1: (color online) (a) Schematic of an optical DFB struc- 
ture consisting of a uniform index grating with a homogeneous 
gain region that realizes a "PT-non-invariant Dirac Hamilto- 
nian. (b) Zeros of M22 (squares) and of Mn (circles) in the 
complex energy plane e = EL for moL — 1 and for increas- 
ing values of normalized gain A. The critical value A c , above 
which bound states emerge and the spectrum of H ceases to 
be real- valued, is A c ~ 1.755/L. (c) Spectral transmission 
\t\ versus normalized energy e of incident wave for increasing 
values of XL: Curve 1, XL = 0; curve 2, XL = 1; curve 3, 
XL = 1.5; curve 4, XL = 1.65; curve 5, XL = 1.74. The 
dotted curve in (c) is the spectral transmission of the DFB 
structure for XL = 1.74 when the gain region is replaced by a 
lossy region. 
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FIG. 2: (color online) (a) Schematic of an optical DFB struc- 
ture consisting of a uniform index grating with two homo- 
geneous gain and lossy regions that realizes a PT-invariant 
Dirac Hamiltonian. (b) Zeros of M22 in the complex energy 
plane e = EL for moL = 1 and for increasing values of A 
(squares: A = 0; triangles: XL = 3; circles: XL = 4.46; 
rhombs: XL — 6). VT symmetry breaking is reached at 
A c ~ 4.46/L. (c) Behavior of spectral transmission \t\ for 
increasing values of AL: Curve 1, XL — 0; curve 2, XL = 3; 
curve 3, XL = 4; curve 4, XL = 4.45. 



side). An example of such a spectral singularity, arising 
from the crossing of an antiresonance with the real en- 
ergy axis, is discussed in the first example below. 

Optical realizations of spectral singularities and VT 
symmetry breaking. Let us specialize the general results 
of the spectral theory by considering two examples of 
complex potentials, in which spectral singularities and 



VT symmetry breaking of the Dirac equation correspond 
to well-known physical phenomena in the theory of DFB 
optical systems. 

As a first example, let us consider a 'PT-non-invariant 
potential corresponding to 9{x) — and m(x) — mo, 
j(x) — —A in the interval |x| < L/2 [Fig. 1(a)], where A is 
the dimensionless gain coefficient. This case corresponds 
to the simplest version of a DFB laser with a uniform in- 
dex grating and a uniform gain region pl| . The transfer 
matrix of this structure is given by 19] 



cosh(pL) — i^ sinh(pL) 



-i r -f- smh(pL) 
cosh(pL) + i - sinh(pL) 



(7) 

where a = —E + iX and p = (toq — a 2 ) 1 / 2 . Note that 
the functional dependence of M on E and A is solely 
via a = -E + iX, so that Mik{E, A) = Mik(E - iX, 0). 
In the Hermitian limit A = 0, the spectrum of H is 
purely continuous, and a number of resonances (i.e. 
zeros of M.22 in the lm(E) < half plane) as well as 
of anti-resonances (i.e. zeros of M.\\ in the lm(E) > 
half plane) do exist [Fig. 1(b)]. As A is increased, the 
spectrum remains real-valued, the resonances (antireso- 
nances) rigidly shift parallel to the imaginary axis, until 
at a critical value A = X c two resonances cross the real 
energy axis, i.e. they become two spectral singularities 
[Fig. 1(b)]. The critical value A c can be expressed in 
the form A c = f(moL)/L, where the function f(moL) 
can be calculated numerically. For instance, for the 
case of Fig.l {m^L = 1) one has X C L ~ 1.755. At 
A > A c the resonances cross the real energy axis and 
bound states with Im(_E) > emerge; correspondingly, 
the spectrum of H ceases to be real-valued. Such a 
transition, from A < A c to A > A c , is accompanied by a 
narrowing of the resonance widths in the transmission 
spectrum as A — > X~ [see Fig. 1(c)] and to the threshold 
for self-oscillation at A = A c . Note that the imaginary 
parts of the eigenvalues E for A > A c are precisely the 
growth rates of the two detuned unstable modes at the 
onset of lasing found in the theory of DFB lasers with 
a uniform grating [l8|, It is remarkable that the 

well-known physical phenomenon of self-oscillation in 
DFB lasers mimics the onset of a spectral singularity 
of a non-Hermitian relativistic wave equation. Let us 
consider now the same DFB structure of Fig. 1(a) but 
with j(x) = A, corresponding to a uniform grating with 
a homogeneous lossy region. The transfer matrix of the 
structure is given again by Eq.(7), but with a = —E — iX. 
In this case, as the loss coefficient A is increased from 
zero, the spectrum remains real-valued shifting parallel 
to the imaginary axis. At the critical value A = A c , 
two antiresonances (rather than resonances) now cross 
the real energy axis, i.e. they become two spectral 
singularities. As opposed to the previous case, at A = A c 
the spectral transmission does not diverge because 
the spectral singularities are born from antiresonances 
(rather than from resonances); see the dotted curve in 
Fig. 1(c). Moreover, the crossing does not correspond, 
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as in the previous case, to the onset of lasing in the 
DFB structure, because at A > A c the two bound states 
supported by the Hamiltonian have now a negative 
growth rate (Im(_B) < 0), i.e. any initial perturbation 
(starting from spontaneous emission noise) is damped 
rather than amplified. 

As a second example, let us propose a "PT-invariant 
DFB structure which could be used to observe VT 
symmetry breaking of the Dirac equation. The struc- 
ture, schematically shown in Fig. 2(a), is composed by a 
uniform index grating with two symmetric homogeneous 
gain and lossy regions, i.e. 9(x) = 0, m(x) — tuq for 
|x| < L/2, j(x) = -A for -L/2 < x < 0, and j(x) = A 
for < x < L/2. The transfer matrix of the structure 
can be calculated as M = M2M1, where Aii and 
M.2 are the transfer matrices of the uniform gain and 
lossy sections [see Eq.(7)]. Figure 2(b) shows the loci of 
zeros of M.22 in the complex energy plane for increasing 
values of XL and for m^L = 1. Note that, owing to 
the PT-invariance of H, Mu(E) = M*2 2 {E*), and 
thus the zeros of M.\\ (not shown in the figure) are 
simply the complex conjugates of those of M.22- As 



A is increased above the critical value A c ~ 4.46/i, a 
VT symmetry breaking occurs, with the appearance of 
two pairs of complex-conjugate eigenvalues belonging to 
the point spectrum of H and arising from two couples 
of resonances and antiresonances crossing the real 
energy axis. As the point of VT symmetry breaking is 
approached, narrowing of the transmission resonances is 
observed [see Fig. 2(c)], and the onset of lasing at A = A+ 
corresponds to the breaking of the VT symmetry. 
Conclusions. DFB optical structures provide a fertile 
classical simulator of non-Hermitian relativistic wave 
equations. This work suggests that well-known phenom- 
ena occurring in DFB structures, like spectral narrowing 
of resonances and self-oscillation, are the measurable 
quantities associated to the onset of spectral singularities 
and VT symmetry breaking of Dirac Hamiltonians with 
certain complex couplings. 
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